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Abstract. We carried out a detailed analysis of spectra of 
the unique roAp star HD 101065 (Przybylski's star) near 
the resonance doublet Li I 6708 A, using a most complete 
line list including all possible transitions between REE 
levels of NIST database. Our model calculations were per- 
formed under two assumptions: blend of Li and REE lines, 
and blend of REE lines only. They prove that Li lines are 
present in the range 6707.72 - 6708.02 A and that the 
resulting Li abundance is 3.1 dex (in the scale log N(H) 
= 12.0), while the isotopic ratio 6 Li/ 7 Li is near to 0.3. 
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1. Introduction 

Among the rapidly oscillating Ap (roAp hereafter) stars, 
HD 101065 is the most unusual (Przybylski HMD EH- 
The report of Przybylski in 1976 at the IAU Symposium 
No 32, showing the very specific properties of the spectrum 
of HD 101065, made a very strong impression on special- 
ists of Ap stars and at once attracted much attention to 
this star. Significant difference between the spectrum of 
this star and the spectra of other Ap stars were pointed 
out by Przybylski : first, the absence or extreme weak- 
ness of lines of elements that appear generally normal or 
close to normal in other Ap stars, like iron and some other 
iron group elements, then the abnormal strengthening of 
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Rare Earth Elements (REE hereafter) lines. Przybylski 
predicted the presence of a magnetic field for this star. He 
proposed to explain the observed anomalies of elemental 
abundances by the stability of the upper atmosphere due 
a strong magnetic field, which would favour radiative dif- 
fusion. Wolff and Hagen (1976) did find a magnetic field of 
about 2.5 KG in this mysterious star. Kurtz and Wegner 
1)1979[1 began photometric observations of HD 101065 in 
the seventies and detected well-defined pulsations with a 
peak-to-peak amplitude of about 0.012 mag and a period 
of 12.14 minutes. 

It was the first Ap star in which pulsations were found. 
Now this phenomenon is known to be typical of other Ap 
stars, which form the group of roAp stars. The first at- 
tempts to identify the lines in the spectrum of HD 101065 
were made by Przybylski Ifl 9B 11 ITT) 6 6|) . and Warner (1966). 
These papers mention for the first time the possible pres- 
ence of the resonance lithium doublet. In the subsequent 
study of this star, Wegner and Petford l|1974f) have pro- 
posed that it is the coolest of Ap stars. They considered 
the very strange spectrum to be the combined result of 
low temperature and high surface REE abundances (the 
abundance excess is about 4 dex). Beginning with Przy- 
bylski, investigators suggested that heavy line blanketing- 
changes the atmosphere structure very strongly, and that 
the computation of a realistic model atmosphere for this 
star is complicated by the necessity to take REE lines 
blanketing into account. Cowley at al. Q1998I I2000fl ex- 
amined in detail the problem of effective temperature to 
define the model atmosphere of HD 101065, and have car- 
ried out a thorough quantitative analysis (abundances of 
54 elements) of the spectrum in the range 3959-6652 A. 
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Cowley and Mathys 1|1998[) have noted in their conclu- 
sion that the presence of the third spectrum REE (Pr III, 
Nd III, Ce III...) is evidence for the unusual structure of 
the atmosphere of this star, which includes a very thin 
surface convective zone typical of stars of spectral type 
F2. One possible explanation of these extreme abundance 
anomalies is radiative diffusion, allowed by the stabiliza- 
tion of the atmosphere by a strong magnetic field. Cowley 
et al. (|19981 I2000f) did not address the question of the 
lithium abundance in the atmosphere of this star. The 
present work is a further investigation of the spectrum of 
HD 101065 in the region 6675-6735 A; the main goal is 
the analysis of the lithium blend 6708 A and estimation 
of the atmospheric lithium abundance. 

The very complicated spectral feature at A 6708 A in 
the observed spectra raises the problem of correct line 
identification, especially REE lines, in this region. A com- 
prehensive analysis of REE lines with latest atomic data 
was performed, and we show that the resonance doublet 
of lithium is the main contributor to the 6708 A feature. 



2. Lithium in the spectrum HD 101065 

The presence of lithium in the spectrum of this unique star 
was first mentioned by Przybylski ( 1961 ), and this remark 
was noticed by Warner (1966}, who made additional ob- 
servations of this star in Radcliffe Observatory, using the 
74" reflector, in the spectral range 3770 - 6880 A (disper- 
sion 6 A/ mm) . He found very strong lines of singly ionized 
REE in the region 5500 - 6880 A, the relative intensities of 
which are similar to laboratory intensities from the tables 
of Meggers et al. {126JJ. However, as mentioned by Przy- 
bylski {MSI), the line Sm II 6707.45 A can be considered 
only as part of the lithium blend, and the main contribu- 
tion to this blend remains the resonant lithium doublet Li 
I 6708 A. The first estimation of the lithium abundance 
relative to the solar value was made by Warner (Q966): 
[Li] = 2.4 dex. This author also mentioned the probable 
presence of 6 Li. 



3. Observations and reduction 

The observations were made by one of us (PN) with the 
1.4m CAT at ESO La Silla in March 1996, in the frame- 
work of the International Project "Lithium in magnetic 
CP stars" . The Coude Echelle Spectrograph was used with 
resolving power R = 100 000 and S/N > 100. The usual 
data reduction (bias and dark subtraction, division by flat- 
field, extraction and wavelength calibration) was done by 
PN during the observing run, using the old IHAP system 
of ESO (see North et al. lT§Ml for more details). In Table 1 
we give the epoch, exposure time and wavelength coverage 
of our observations. 



4. Model Atmosphere 

We have computed a model atmosphere with T ff — 
6600 K, log g = 4.2 and abundances of most (54) elements 
from Cowley et al. I|2000J) . The other element abundances 
were adopted from Kurucz lfHT9 9). A detailed discussion 
of model atmosphere parameters for HD 101065 can be 
found in the above mentioned paper, as well as in the 
paper by Piskunov and Kupka (2001). 

The most recent method for determining fundamental 
stellar parameters through pulsation analysis was applied 
by Matthews et al. (0,999) to some roAp stars, including 
HD 101065. But it needs a photometric calibration. Hp 
calibration results in T e s = 7400 K, which is too high, as 
confirmed by detailed spectral analysis (see Cowley et al. 
|2TJUU|). We used the SAM 12 program (see Pavlenko l20CJ3al 
I2003bjl . a modified version of Kurucz's code ATLAS12 
(Kurucz 119991) . The essential changes were made in the 
subroutines accounting for opacities. Namely, the opacity 
caused by absorption of atomic lines was taken into ac- 
count with the Opacity Sampling (OS) technique (Sncden 
et al. I1976|) . Then, bound- free absorption cross-section of 
C I, N I and O I atoms in the wide spectral region 30- 
8880 nm were calculated by us, Pavlenko & Zhukovskaya 
127)011 using the data of the TOPBASE project (Seaton 
1992). New opacity tables of the bound- free ab- 
sorption of C I, N I, O I atoms, reduced to the 
format of Hofsaess (1979) opacity tables, are avail- 
able on the web site of Ya.P.( 2003c). The importance of bf 
absorption for formation of spectra of different stars are 
discussed elsewhere (PavlenkoiTffil^l l2003al Asplund et al. 
ISUUUjl . 

The basic advantage of the SAM12 code is the possibil- 
ity to calculate model atmospheres for stars with chemical 
composition different from the solar one, as well as in AT- 
LAS^. Abundances of elements were set as input data 
which do not change with depth. The microturbulence ve- 
locity was set constant with depth. 

The profile of the line absorption coefficient was 
defined by a Voigt function H(a,v); the damping 
constants were taken from line databases or calcu- 
lated within the frame of the Unsold approximation 
(Unsold 1955). The list of absorption lines of atoms 
and ions for the range AA40 — 60000 nm was compiled 
using the VALD database (Kupka et al. 1999); the 
REE lines data, calculated with Cowan's code (Quinet 
et al. I1999|l . are taken from the DREAM database 
(http://www.umh.ac.be/' astro / dream, shtml \. Unfortu- 



nately, this database contains lines of only about one 
quarter of all rare earth ions. Therefore, in the opacity 
sampling calculation we consider two samples of REE 
lines. In the first we calculated the absorption provided 
by the VALD and DREAM lines. Then, we amplified the 
absorption of DREAM lines by a factor of 10 (in addition 
to the VALD lines) to compensate for the incompleteness 
of the REE line list. The comparison of these two cases 
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Table 1. List of observations for HD 101065. 





Date 


UT 


Exp. 


HJD 


Range 


# 


d.m.y 


h m 


[m] 


2450000+ 




04 


11.3.96 


5 11 


20 


153.726 


6675-6735 A 


26 


10.3.96 


4 41 


20 


152.705 


6120-6180 A 



shows that REE lines opacity does not affect significantly 
the temperature distribution in our model atmospheres. 

Kurucz's models with T cS = 6750 K, logg = 4.0 (Ku- 
rucz I1999f) and abundances of Cowley et al. I|2000|l were 
also used in calculations of synthetic spectra. We calcu- 
lated 6750/4.0/[M] = 6750/4.0/0 model with solar abun- 
dances, as described previously, and compared them with 
Kurucz models (Kurucz 1993). Comparison of tempera- 
ture structures and corresponding synthetic spectra shows 
their similarity. As far as our model atmosphere is rather 
iron deficient, the increase of opacity due to the REE ele- 
ments may be compensated by the weakness of iron lines. 

Cowley et al. analysed the spectrum of HD 101065 
with ATLAS-9 model with Opacity Distribution Function 
(ODF) tables, calculated by the authors on the basis of 
scaled solar abundances to fit better the observed color in- 
dices, in particular the unusually low Ci index. The model 
atmosphere for this star was also computed by Piskunov 
and Kupka l|2001j> . also through ODF method, with a 2.25 
dex Fe overabundance and a 2.49 dex Ca excess to simu- 
late unknown opacity in ODF tables. These authors used 
the abundances of Cowley et al.JSUUUJ) for model computa- 
tion with these ODF. We have followed this approach for 
OS method, by computing line opacity due Fe I and Ca I 
lines with factors for gf- values of 177.8 (2.25 dex) and 308 
(2.49 dex) respectively. In Fig. 1 we compare the T-P rela- 
tion for model atmospheres with T c ff = 6750 and log g = 
4.0 of Kurucz (solar abundance and [Al] = +1) and ours 
with Cowley's abundances, calculated in two ways - one 
OS model, similar to the ODF model of Piskunov & Kupka 
(SAM12:6750/4.0/0, analogue PK2000), and our 
two OS model atmospheres for T c s = 6750 and 6600 K 
and log g — 4.0, with gf(REE) lines enhanced by a fac- 
tor of 10(SAM12:6750/4.0/0/Cowley 2000, g/ ree *10 
and SAM12:6600/4.2/0/Cowley 2000, 5 / ree 2*10). 
In general, the photospheres of models with higher opaci- 
ties are shifted upwards, into the layers of lower pressure; 
however, the differences in the T-P relation are rather 
small. 

It is worth noting, that our OS model atmospheres 
6750/4.0/0 and 6750/4. 0/+1 agree with corresponding 
Kurucz's ODF models within 50 K in the line forming re- 
gions. Therefore, our OS model atmosphere with "metal 
rich" opacities should follow the Piskunov & Kupka (2001) 
ODF model. 

To analyze the Li blend, we used an atmosphere model 
with T c ff = 6600, logy = 4.2. The resulting abundances 



were compared for this model and for the OS model, which 
is similar to that of Piskunov & Kupka for the same pa- 
rameters, both for Cowley's abundances, and for Kurucz's 
model for solar abundances. They are cited and discussed 
in the next section. 



5. Li blend analysis and REE lines contribution 
to the 6705.75-6708.75 A spectral region 

We calculated the synthetic spectra in the range 6705.75- 
6708.75 A with the help of the STARS P ( Ts ymbal lTTO 
and ROTATE codes (see Shavrina et al. 200(|, both devel- 
oped by V. Tsymbal, and of the URAN code (Yushchenko 
1997|. All cited codes are based on Kurucz' modules for 
synthetic spectra ( " S YNTHE" , GHH . The URAN pro- 
gram allows to select the elemental abundances automat- 
ically, and applying it to the available spectral ranges 
(6120-6180, 6675-6735 A), one obtains results which are 
close to those of Cowley et al.(2000). 

The line list for spectra 
was compiled from DREAM 



calculations 
database 

{http://www.umh.ac.be/' astro / 'dream. shtml\ and VALD 
CKupka et al. I1999fl . When there was a choice, the 
preference was given to gf values from DREAM. The 
wavelength and gf-values for lithium lines ( 7 Li and 6 Li) 
were taken from the paper of Smith et al. 0998). In addi- 
tion, we have calculated wavelengths of REE lines, singly 
and doubly ionised, using the data of REE energy levels 
from the NIST database l|2001|l according to selection 
rules ( 8 J = 0, ± 1) and combining even and odd levels 
of energy (odd-even or even-odd). All possible transition 
between REE energy levels of NIST were included in the 
line list for synthetic spectra calculations. For additional 
to VALD and DREAM lines which are rather distant 
from Li lines and do not contribute noticeably to the 
6708 A blend ( see Table 4), ^/-values were matched for 
better agreement of observed and calculated spectra. 

g /-values for the two lines nearest to the lithium lines 
(Tables 2, 4), Nd II 6707.755 A and Nd II 6708.030 A, 
were estimated by P. Quinct with Cowan's code using for 
the second line the trial identification of the top level, 
which was unspecified in the NIST list of REE energy lev- 
els. This line, Nd II 6708.030 A forms a blend with the 
strong line Ce II 6708.099 A, whose gf value is contained 
in the DREAM database. In the region studied, there is a 
rather pure line Ce II 6706.051 A . DREAM value of logg/ 
= -0.95 for this line provides Ce abundance log N(Ce)= 
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Fig. 1. T-P relation for model atmospheres considered in this work. Two models from the Kurucz's (1993) 
grid are showed, as well as SAM12 simulation of model atmosphere using Piskunov & Kupka (2000) 
approach; SAM12 model atmosphere computed with Cowley (2000) abundances and REE bf opacities 
enhanced by factor 10 of the same T ffand log g; structure of model atmosphere T e s= 6600/4.2 with 
Cowley abundances and enhanced REE opacities is shown also. See text for more details. 



-7.63, which almost coincides with the estimate of Cow- 
ley et al. (HHP, log N(Ce) = -7.60. This Ce abundance 
value results in estimate log N(Nd) = -7.97 from the line 
Nd II 6708.030 A with log 5/ = -1.13. This estimate dif- 
fers from the value of Cowley et al. 1)2000(1 log N = -7.65 
(±0.28) by 0.32 dex, but it only slightly exceeds the er- 
ror value. Therefore, Quinet's estimate for the gf value 
of Nd II 6708.030 seems a reasonable one. However, its 
blending with the 6 Li lines leads to some uncertainty in 
the estimation of the lithium isotopic ratio. The line of 
Sm II, A 6707.779, almost coincides with the centre of the 
lithium blend, but unfortunately, the calculations for Sm 
II are not possible because the energy matrix dimensions 
exceed the permitted sizes for Cowan's code. We tried to 
calculate the blend 6708 A with Sm II line absorption A 
6707.779, instead of Li lines, adjusting the gf value for 
better agreement with the observed blend. In this way 
we could estimate an upper limit to log gf for this line, 



—0.63. The resulting fit was bad (see Fig. 2). Only the 
inclusion of 7 Li and 6 Li lines (see Table 2) with the ratio 
6 Lij 1 Li = 0.3 allowed us to represent the observed profile 
of blend rather well (see Fig. 2). Our estimate of log gf for 
Sm II line A 6707.779 in this case is equal to -2.68. We 
should repeat, that our estimate of lithium isotopic ratio 
is influenced by blending with the Nd II line 6708.030 A. 

Table 2 contains the line list used in model calcula- 
tions. The wavelengths of lines marked by asterisk were 
calculated on the basis of NIST energy levels (2001). Cor- 
responding g /-values were estimated by us using element 
abundances from Cowley et al. ( 20QHJ), corrected for some 
REE elements (see Table 3) on the basis of lines with 
known gf- values (from DREAM and VALD lists). We 
used a value of microturbulent velocity V t = 2 km/s to 
compensate for magnetic line broadening. This value was 
selected by us for better agreement of computed and ob- 
served line profiles. 
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Table 2. The list of lines used for spectra calculations in the range 6705.75 - 6708.75 A 



el 


A,A 


E" 


logs/ 


source 




el 


A,A 


E" 


log gf 


source 


*Dy II 


6705.727 


2.078 


-2.68 






* Gd II 


6707.462 


3.270 


-1.98 




Nd II 


6705.891 


3.269 


-2.70 






Sm II 


6707.473 


0.930 


-1.48 


VALD 


Ybll 


6705.965 


5.856 


-3.04 


DREAM 




Yb II 


6707.603 


6.651 


-1.38 


DREAM 


Ce II 


6706.051 


1.840 


-0.95 


DREAM 




* Sm II 


6707.648 


1.746 


-1.27 




Tm II 


6706.150 


5.322 


-1.07 


DREAM 




Nd II 


6707.755 


0.170 


-3.55 


DREAM [Quinet] 


Tm II 


6706.262 


3.955 


-2.36 


DREAM 




7 LiI 


6707.756 


0.000 


-0.427 


[Smith et al.] 


Ce II 


6706.307 


3.195 


-2.40 


DREAM 




7 LiI 


6707.768 


0.000 


-0.206 




Pr III 


6706.492 


3.104 


-1.28 






*Sm II 


6707.779 


2.037 


-2.68 




Pr III 


6706.705 


0.550 


-1.64 


DREAM 




7 LiI 


6707.907 


0.000 


-0.932 


[Smith et al.] 


Nd II 


6706.738 


2.868 


-2.48 






7 LiI 


6707.908 


0.000 


-1.161 


55 


Sm II 


6706.789 


1.586 


-2.00 






7 LiI 


6707.919 


0.000 


-0.712 


!) 


Sm II 


6706.807 


1.874 


-1.78 






6 LiI 


6707.920 


0.000 


-0.478 


55 


Tm II 


6706.906 


4.908 


-2.47 


DREAM 




7 LiI 


6707.920 


0.000 


-0.931 


!; 


Nd II 


6706.922 


3.211 


-0.88 






6 LiI 


6707.923 


0.000 


-0.179 




Nd II 


6707.015 


1.490 


-1.88 






Nd II 


6708.030 


1.522 


-1.13 


DREAM [Quinet] 


Nd II 


6707.033 


2.222 


-3.68 






6 Li I 


6708.073 


0.000 


-0.304 


[Smith et al.] 


Ce II 


6707.121 


1.255 


-3.76 


DREAM 




Ce II 


6708.077 


2.250 


-2.57 


DREAM 


Dy II 


6707.153 


3.292 


-1.27 






*Er II 


6708.088 


3.155 


-2.58 




Dy II 


6707.266 


2.890 


-1.28 






Ce II 


6708.099 


0.701 


-2.12 


DREAM 


Sm II 


6707.342 


0.884 


-2.00 






*Nd II 


6708.400 


3.192 


-2.48 




Er II 


6707.418 


3.482 


-1.44 






*Nd II 


6708.458 


3.536 


-1.08 




Nd II 


6707.433 


1.499 


-2.17 






*Nd II 


6708.629 


0.746 


-4.58 




Nd II 


6707.453 


2.880 


-3.18 















*g/-values were estimated by us from best fitting with observed spectrum. 



In Table 3 we present our results for Kurucz's and our 
models with gf(REE)*10 and OS model (PK in the table) 
which is similar to that of Piskunov and Kupka (2001), 
and compare them with data of Cowley et al.( 2000) for 
Vt = 2 km/s. A better agreement was achieved for Ce II on 
the basis of several Ce II lines with DREAM gf- values, the 
lithium abundances being almost the same for all models 
used in this work. We were forced to reduce the Nd II and 
Sm II abundances taken from Cowley et al. ( 2000) when 
we fitted the observed line profiles with theoretical ones. 
These values (see Table 3), based only on one line of each 
element, 6708.030 A of Nd II and 6707.473 A of Sm II, 
can not be considered as reliable abundance estimates. 

We should note that the distinctive feature of this star 
among other roAp stars is the predominance of numer- 
ous strong Ce II lines in the spectrum, which betrays the 
extreme overabundance of this ion. Not accounting for 
DREAM line list of Ce II can result in overabundances 
for some REE and other elements. 

In Table 4 we show the shares of the main contributors 
in total absorption (in percent) in each wavelength with a 
0.02 A step, as they are calculated by the STARSP code, 
for the lithium doublet range 6707.60 - 6708.16 A with- 
out instrumental smoothing, for the model atmosphere 
6600/4.2 with V t = 2^8-! and 6 Li/ 7 Li = 0.3. 

Reyniers et al. ( 2002 ) draw attention to a possible error 
in Li line identification (instead of Ce II line 6708.03) in 
the spectra of s-process enriched post-AGB stars. In our 
high resolution spectra, these two lines are sufficiently far 



apart to be distinguished, as can be seen in Figure 2 and 
Table 4. 

6. Discussion of the results 

Since the discovery of the star HD 101065 by Przybyl- 
ski, the debates relating to the effective temperature of 
this star have not settled yet. Apparently, the reason lies 
in the large blanketing effect due to strong absorption in 
REE lines, identified and unidentified. Incompleteness of 
data for these lines does not allow to take into account cor- 
rectly this effect in model atmosphere calculations. Cow- 
ley et al. ( 2000 ) have clearly shown a large excess of REE 
abundances ( 4 dex ) and deficiency of iron peak elements 
( Fe, Ni ), while Co has appeared in excess of 1.5 dex. The 
list of REE lines in this work included new identified lines 
in the spectral range bf 3959-6652 A. Our results include 
the range 6670 - 6735 A where lies the remarkable feature 
with the resonant lithium doublet line 6708 A. 

We have carried out detailed calculations of the blend 
6708 A(6705.75 - 6708.75 A), using the REE list of 
atomic data from the DREAM and VALD databases. 
The evidence for presence of lithium in the spectrum of 
HD 101065 is most probable, as shown by the excellent fit 
of synthetic spectra which include the Li lines, while the 
synthetic spectra without lithium (only REE lines, Sm II 
6707.799) clearly fail to achieve a good match. All possi- 
ble transitions between REE energy levels of NIST were 
included in the line list for synthetic spectra calculations. 
The presence of other REE lines near 6708 A, that would 
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Table 3. Element abundances for modelling 6705.75-6708.75 A spectrum 



el 


'Li only 


b Li/'Li = 0.33 






6600/4.2 


6600/4.2 


6750/4.0 


6600/4.2(PK) 


6600/4.2(Cowley et al.( 20001 


Li I 


-8.90 


-8.90 


-8.80 


-8.83 




Ce II 


-7.63 


-7.63 


-7.65 


-7.63 


-7.60+-0.26 


Pr III 


-7.55 


-7.55 


-7.60 


-7.65 


-7.46+-0.16 


Nd II 


-7.97 


-8.07 


-7.87 


-8.07 


-7.65+-0.28 


Sm II 


-8.32 


-8.32 


-8.25 


-8.42 


-7.75+-0.19 



Table 4. The shares of main absorption contributors in the range 6707.60 - 6708.16 A 
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not be included in our list remains improbable. All tran- 
sitions near 6708 A are due low energy levels, and low 
energy REE levels that would not belong to NIST are un- 
likely. The only line which might mimick the Li lines is Sm 
II 6707.799 A . But we showed that absorption in this line 
is not large enough to account for the observed spectral 
feature, while a perfect fit is obtained when accounting for 
the Li line components, which are distributed on a rather 
wide range (see Tables 2 and 4). 

In our work the choice of parameters of synthetic spec- 
trum near the line Li I 6708 A is described in detail taking 
into account the REE lines with known and unknown gf 
values. The variants of the calculations, with and without 
lithium lines, are considered, which prove the presence of 
a lithium doublet 6708 A in the spectrum of HD 101065 
with a lithium abundance of 3.1 dex ( in the scale log 
N(H) = 12.0 ). The observed blend is best represented 
at the isotopic ratio 6 Li/ 7 Li = 0.3. The latter value can 
be considered as an upper limit in the current state of 
atomic data; more accurate estimation of the lithium iso- 
topic ratio would require better spectral resolution and 
better ^/-values for REE lines. 

We carried out the calculations with Zeeman and 
Paschen-Back splitting of the Li line. For H=2300 G, the 
Paschen-Back effect is negligible (this is common for mag- 
netic fields lower than 4000 G - see Shavrina et al. I2000J) 
and can not induce any additional asymmetry. Zeeman 



broadening is also unimportant, due to the significant sep- 
aration (0.15 A) of the components of the Li doublet. On 
the basis of our spectra calculations in two wider spec- 
tral ranges, 6120 - 6180 A and 6675 - 6735 A, the value 
V t = 2 km/s was adopted, which must compensate for the 
magnetic broadening of other lines. We have also consid- 
ered the possible binarity of HD 101065. This problem was 
discussed by Hubrig et al. (12000|l. They noticed that none 
of the 14 roAp considered, including HD 101065, is known 
to be a spectroscopic binary (although several are mem- 
bers of wide visual binary systems); therefore the 6708 A 
feature cannot be due to a secondary component of this 
star . 



Warner's ( 1966) estimate of the lithium abundance for 
this star from a spectrum with 0.3 A resolution was 2.4 dex 
relative to the solar value, that is, 3.3 dex relative to hy- 
drogen. He also assumed the presence of 6 Li on the sur- 
face of the star. The good agreement between our results 
and those of Warner (19@ is encouraging. Warner had 
tried to explain the large excess of rare earth elements 
and the peculiar structure of the atmosphere by dragging 
out r-process elements from deeper to surface layers of 
the star. However, the high lithium content observed, 3.1 
dex ( on the scale where log N(H) = 12 ), is not unex- 
pected for normal F stars, which are separated on high 
and low Li abundances (Chen et al.. l20()l)l . while Warner's 
assumption would imply complete mixing and, therefore, 
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Fig. 2. The fitting of observed and calculated spectra of HD 101065 near 6708 A : dots: observed spectrum; dashed 
line: calculated spectrum taking into account lines of the main isotope 7 Li only; thick line: spectrum with the ratio 
e Li/ 7 Li = 0.3. The thin line corresponds to a spectrum computed without Li but with the Sm II line 6707.799 A. The 
positions of those lines which are the main contributors in absorption, are marked at the top of the figure. 



complete lithium destruction. Thus the radiative diffusion 
mechanism seems a much more promising explanation, es- 
pecially as the magnetic field, reaching 2300 G (Cowley et 
al.. I2000|l . would make it easier. On the other hand, the 
isotopic ratio 6 Li/ 7 Li = 0.3 might be considered as indi- 
cating spallation reactions on the stellar surface, though 
a combination of radiative and ambipolar diffusion might 
perhaps lead to the same result. 

The Hipparcos measurements have confirmed that the 
star is a dwarf or a subgiant: using its parallax (7.95 ± 
1.07 mas) and the adopted T c g = 6750 K, we obtain 
logg = 4.00±0.13 by interpolation in evolutionary tracks, 



assuming a normal evolution and a solar global composi- 
tion. However, its spectrum is reminiscent of S-stars (ex- 
cept for Li !), whose atmospheres are enriched by fresh 
synthesized material from deep layers. Whether the ex- 
treme anomaly of the spectrum of this interesting star is 
just another case of radiative diffusion, or requires a com- 
pletely new explanation, remains to be seen. 
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